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Abstract 
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I. INTRODUCTION 



Quite a number of works were recently devoted to the studies of optical properties of 
nanodimensional bodies. Special attention has been paid to the metallic nanoparticles, with 
the help of which one could amplify the near electric field at the frequencies of the localized 
plasmon resonances jl| • Possible applications of that effect have been considered. The most 
developed of them is the use of large local fields near a rough surface for increasing the 
surface-enhanced Raman scattering (SERS) Modification of properties of the radiating 
atoms by nanobodies of different shape and composition has been the basis for research and 
development of nanobiosensors ^-[sj], nanolasers Q], microscopes that are ca pab le to detect 



single molecules 
others. 



, the devices for DNA decoding [9|, chemical sensors [101 ]- 



121 ]. and many 



Very recently, the plasmon nanostructures were successfully used in optical devices, in 
which the diffraction limit had been overcome thus providing an effective interaction between 

iQ.—nnasH-H. Studying 



near and far fields, e.g. in hyper- and nanolenses 
optical properties of the plasmon nanoparticles made up a new line of investigation such as 
"plasmonics" or "nanoplasmonics" 21 



The above-mentioned applications are based on properties of single nanoparticles or clus 



ters of such nanoparticles. Today, the optical properties of metallic nanospheres 
nanospheroids 3l]- 34|, and some other nanobodies 



21 



22 



23|-|30|, 



351 ] . have been sufficiently well 



studied with the help of known analytical approaches. Many of the other shapes of nanopar- 
ticles have been studied only numerically, which is often insufficient for understanding of 
underlying physics and for development of new nanodevices. 

It is astonishing but an important class of nontrivial nanoparticle shapes , which can be 
described analytically, turns out to be almost uninvestigated. We keep in mind nanoparticles 



of triaxial ellipsoid shape. Such nanoparticles have very different shapes and are o 



ied experimentally 36]- 39|. Despite a great number of theoretical works 40]- 52] devoted 



ten stud- 



to that problem, there had not yet been achieved any appreciable progress in studies of 
optical properties of triaxial ellipsoids. This is due, first of all, to an extreme mathematical 
complexity of the problem. In practice, to describe optical properties of triaxial ellipsoid 
nano par ticles one makes use of a simple solution for a nanoellipsoid in a uniform electric 
field [53] that describes only the simplest ( "dipole" ) plasmon oscillations. Such a description 



2 



is sufficient if scattering of plane electromagnetic waves by nanoparticle is of interest. But 
for many nanooptical problems there are no plane waves and such a solution cannot be used. 
In particular, it does not work when considering an ellipsoidal nanoparticle interaction with 
highly nonuniform fields, which is fundamentally important for a nanoworld. For exam- 
ple, the radiation fields of atoms, molecules, or quantum dots placed near the nanobodies 
are essentially not uniform, and they effectively excite not only "dipole", but more com- 
plex plasmon oscillations also. For example, exactly these modes explain very important 
phenomenon of radiation quenching. 

The aim of this work is to study optical properties of a triaxial metallic nanoellipsoid 
within the quasistatic approximation, in respect to the excitation of high-order plasmon 
oscillations that are somehow analogous to the quadrupole, octopole, and hexadecapole 
spherical harmonics. Here we obtain explicit analytical expressions, which is especially 
important for fast calculations and estimations. 

The content of the remaining part of the article is as follows. In section II, general 
approach to finding the eigen-functions of plasmon oscillation electric field potential and the 
permittivity (frequencies) corresponding to those oscillations is described. In section III, the 
above method is applied to determine explicit values of the permittivity and the plasmon 
mode potentials. In section IV, the graphic illustrations of the found solutions are given and 
the results obtained are discussed. The geometry of the problem in question is illustrated 
in Fig.l. 



II. GENERAL METHOD OF SOLUTION 



The localized plasmons naturally arise as the nontrivial solutions of a homogeneous qua- 
sistatic problem, that is, the Laplace equations with the corresponding boundary conditions 




an 



e— 
On 



(1) 



where e is the permittivity of space around a nanoparticle that we take equal to 1; 
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g n , the plasmon potentials inside and outside a nanoparticle, respectively; -^f L stands 
for a normal derivative of the plasmon potential <p n at a particle surface; e n , the resonant 
permittivity corresponding to the plasmon oscillations of the order n, and A, the Laplace 
operator. The latter equation in ([!]) stands for a fulfillment of continuity of the induction 



normal components . One may show [54], |55j that a homogeneous system ([I]) has nontrivial 
solutions (localized plasmons) for some (negative) parameters e n only. It is naturally to 
call e n as resonant permittivities. If the permittivity of a real particle is close to one of the 
values e n then it will have a plasmon resonance at frequency u n , at which 

(w B ) • (2) 

That frequency uo n is called the plasmon frequency and can be found if the dispersion law 
e (uj) is known. In what follows we will use the Drude's dispersion law, e (uj) = 1 — (u pl / 'uj) 2 , 
where uj v \ is the bulk plasma frequency of nanoparticle material (metal). 

The eigenfunctions of the plasmons e n (r) = — V0 n (r) (where V is the gradient operator) 
are orthogonal in sense that (n 7^ m) 



dr (e n ■ e m ) = (3) 



v 



These eigenfunctions make the complete system, and the solution for a nanoparticle problem 
in external field E° may be expanded over wavefunctions e n (r) of the localized plasmons 



55] 



E (r) = E° (r) + $> (r) ( I^L± ) ^ (4) 

^— ' \e n - e (u) J J dr (e n ■ e n ) 

v 

where e (uj) is the dielectric constant of nanoparticle and the integration is performed over 
a nanoparticle volume. Thus, this problem is reduced to finding the potentials of plasmon 
modes and corresponding eigenvalues (resonant permittivities). 

To determine the properties of localized plasmons in a triaxial nanoellipsoid it is conve- 



nient to use an integral formulation, instead of system ([!]) (see, for instance, 56] ) 



47T 



J dr'((V'^(r'))-V)-^-, (5) 



v 
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where the integration is performed over the nanoparticle volume; r', r are the coordinates 
of some two points inside the nanoparticle, and V, V, the operator gradient over the 
coordinates r', r, respectively. Because the r point is inside ellipsoid, expression (jSJ) is the 
integral equation. The solution of (jSJ) gives us both the eigenvalues e n and the values of the 
plasmon potentials at any spatial point inside ellipsoid. 

Due to eq. ([T]), the solution of integral equation (jSJ) is a bounded harmonic function of 
ellipsoidal symmetry ( ellipsoidal harmonics). The latter may be represented in the form 



of the Niven's functions 

{x 1 ,x 2 ,x 3 y. 



571 ]- 591 ] which are the polynomials of the Cartesian coordinates 



Xi X 2 X 3 

®(x 1 ,X 2 ,X 3 ) = < 1 x 2 X Z Xi X X X 2 X 3 

xGi (x 1 ,x 2 , x 3 ) 6 2 {xi,x 2 , x 3 ) ...O p {xx,x 2 , x 3 ) 



where 



(6) 



Op (xi,x 2 ,x 3 ) 



Xi 



+ 



.1 n 



+ 



6 P a\ + 6 P a 3 + Op— 



- 1. 



(7) 



In Eq.([7j), Op is a constant. By using the ellipsoid coordinates [57j it is easy to prove that 
the Niven's functions are the eigenpotentials of plasmons inside an ellipsoid. 

Thus, knowing the eigenfunctions ([6]) of eqs. ([1]) and (jSJ) one can represent the corre- 
sponding eigenvalues ( resonant values of permittivity) in the following form 



— ^ / dr'((V'0n(r'))-V)-i 
4vr0 n (r) J |r' - 



(8) 



and the problem is thereby reduced to calculating the integral of the Niven's functions ([6 
over a nanoellipsoid volume. The integration can be made by using the Dyson's theorem 60, 
611 ] (see Appendix A). According to it, an integral over an ellipsoid volume of a polynomial 
times Green function 1/ |r' — r|, proves to be the polynomial as well, whose power is greater 
by two 



/ 

v 



dt 



I^X 2 ^3 



Pp+q+s+2 (Xl, X 2 , X3) , 



(9) 
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where P p+g+s+2 (xi, x 2 , x 3 ) is the polynomial of the Cartesian coordinates ( ^ (x a /a a ) < 1) 

a=l 

of the power p + q + s + 2. 

To find a solution for the plasmonic oscillation potential outside an ellipsoid one can use 
again ([S]) 

9n (r) = ^ | dv> ((VV„ (r')) • V) pl^. (10) 

v 

but in eq. (fT0|) . unlike ([5]), r are coordinates of the point outside the ellipsoid. Knowing 
the electric potential of plasmon oscillations inside and outside a nanoellipsoid, as well as 
the corresponding values of the permittivity e n , one can find electric fields near and inside 
ellipsoidal nanoparticles in an arbitrary external field (jl]). 

III. FUNCTIONS OF ELECTRIC POTENTIAL AND VALUES OF PERMITTIV- 
ITY CORRESPONDING TO PLASMON OSCILLATIONS OF A TRIAXIAL NA- 
NOELLIPSOID 

In this section, we will apply our general method to find the plasmonic modes of the 
power n = 1, 2, 3, 4. 

The case n = 1 "dipole" modes 

Here, in accordance with ([6]), the potential inside the ellipsoid is linear over the coordi- 
nates 

1 = ^4 Q) a; Q . (11) 

a=l 

By substituting ( fill into ([5]), and making the required differentiation and integration by the 
Dyson's theorem (see Appendix A) we obtain the following system of equations: 

A[ a) = -^(e 1 -l)a l a 2 a 3 I a A[ a \ (12) 

oo 3 

where I a = J ^ U+ ^ R ^ an d R(u) = Yl ^/ u + a a- As follows from ([121) . there are three 

o Q ct=l 

different eigen-functions (a = 1,2,3) describing the electric potential of the plasmonic os- 
cillations 



6 



= *a, (13) 

and three corresponding values of the permittivity 

e W=l-(^Lj a )- 1 . (14) 

The found eigenfunctions (fT3l) and the eigenvalues ffl4|) are well known and widely used for 
description of the plasmon resonances in ellipsoid nanoparticles in uniform field [53 1. 

The electric potential of the plasmons with n — 1 outside the particle can be found from 
ffTUl) by using ([131) and jUj) (see Appendix B). 

The case n = 2 "quadrupole" modes 

Here, the ellipsoid potential is described by the Niven's functions of two types. Functions 
of the first type are proportional to 0i (x±, X2, x^), that is, can be represented in the form: 



/3=1 Z 3 ,8=1 



where A' and A" are the roots of a quadratic equation tv-a? = ^' ^at * s ' 











HI 











3 / 2 2 2 

^4 a 1 a 2 Q,3 



E 



3 

a=l 

1/2 3 



a 2 



£4 • ( 16 ) 



a=l 



By substituting (fl5|) into (jHJ) and calculating the emerging integrals we obtain explicit ex- 
pressions for the permittivity of those modes: 



f3 2 t 

2 ^-f iV~^a^ 



\ a=l 



/ 3 2 r 

(2) _ 1 . / QlQ2Q3 \ ^ flq-lo 

£ 2 - 1 + ~ / > T7, : 



-1 



»1»2"3 \^ Va /-.„>. 
" I 9 A //_ a 2 I • HH 

a=l c 



Another type of the solutions for the plasmonic modes with n = 2 includes three Niven's 
functions 

# = Z 2 *3, 4 5) =^X 2 . (18) 
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By substituting (1TB]) into (jSJ), and using the Dyson's theorem one can find the corresponding 
values of the dielectric permittivity: 



e <3, = 1 _^ (a 2 + al)/23 )- 

4 4 > = i - {^p (4 + 4) & 

(5) _ 1 _ [010203 ( n 2 



1 



l-PP(a? + «?)M , (19) 



where /„„ = J ,„ + „ S)( ^ )BM («■<?= 1,2,3). 

The plasmonic oscillation potentials outside nanoparticle can again be found by using 
(TTO]) (see Appendix B). 

The case n = 3 "octupole" modes 

In this case, the ellipsoid potential is also described by the Niven's functions of two types. 
Functions of the first type have the form (a = 1, 2, 3) 

^(sA4 (2o) 

where A' and A'' are the roots of the quadratic equation V = 0, that is, 

0=1 a a f3 



Ik) 






HI 











1 6 



4^aJ-3a^-(a 1 a 2 a 3 ) 2 ( ^ ^ + 2/ 

/3=1 \/3=l °^ . 



3 



+2J24-*1 ■ (21) 

/3=1 J 

By substituting (120!) into ([5]), and making the required differentiation and integration with 
the help of the Dyson's theorem one can find the following expressions for the corresponding 
resonant dielectric permittivity 



_(2a-l) 



010203 / A , \ (1 + 2<5 a/9 ) apap 

A \ a a) a / _ 2 

3 



^3 



(2a) 



1+ ^ (A „ + 2a , )E (i + 2M4^ _ (22) 



A'' - o% 
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where 5 a p is the Kronecker's delta-symbol equal to unity at a = j3, and equal to zero at 
[3. 

One more type of the plasmons for n = 3 is due to only one function of the form 

03 ?) = £ix 2 x 3 . (23) 

By substituting (123]) into (jHJ), and using the Dyson's theorem one can find the corresponding 
values of the permittivity 

^ = l-\^^-(j:<^)lm)\ (24) 



V Q = 1 



00 

where J 123 = / 


The plasmon electric potentials outside ellipsoid can be found from (JTUl) by using ( I20l) - (l24l) 
(see Appendix B). 

The case n = 4 "hexadecapole" modes 

Now, there are nine different functions of the plasmon oscillation potentials. Six of them 
are as follows (a = 1,2,3) 



,(2a-i) _ 



3 T 2 

j2 , p 2 + 1 

3 „2 



3 ^ 3_„2(5 

where and Vt'^ are the roots of equation V = 0: 

0=1 "/3 



O" 



1/2 



4{ 4 (E4-^j-9a^ ^4-2^) +^[ . (26) 
The rest of the plasmon eigenf unctions are (a = 1, 2, 3) 
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3 2 \ / 3 2 



v /3=l ^ / V/3=l " a Z 3 / 

where and 2^ are three different pairs of the roots of the system of equations: 



3 1 

Est 1 



a=l 
3 



al E'-E" 



E-A 2 - + ^ ± -7 = °- (28) 



c" - a 



The expressions for 5^, and are rather complicated, so we do not quote them in the 
explicit form. 

The corresponding to (1251) values of the permittivity are represented as 



-(i) 



.(3) 



1 + 



aia 2 a 3 ^ Q[ - a 2 a 

E(l + 25q, 2 + 2<5 a3 ) (a 3 J a2 + all a3 ) , 2 2 x 
— 2 (a 2 + a 3 ) J 23 

4 J / 2 , 2^2 (1 + 2(5 a3 + 2jgi) a 2 J a zx 

\ (a3 l] ^ 



(2i(2 2 Q. 3 
3 



E(l + 25 a3 + 2S al ) a 2 a {ajl a3 + a§ lal) , / > > , 
2(a 3 + a 1 )/ 31 



a=l 



-(5) 



OlQ. 2 (l 3 



{/ 2 , 2^ 2 (1 + ^^Ql + 2^a2) CLgIgl2 

a=l 3 a 



3 ^ _1 

^ (1 + 25 a i + 2<5 q2 ) a* (ap Q i + a? -^a2) „ / 2 , 2\ T \ 
- £l' 3 -a 2 2 (a x + a 2 J i 12 f , (29) 



a=l 



where J a/?7 = / (u+ ^ )(tH , B ^ (ll+ ^ )fl(tl) (a, ft 7 = 1, 2, 3). Expressions for ef\ ef\ ef can be 
obtained via a substitution Q' a — > from the expressions for e 4 , £4 , £ 4 , respectively 

The remaining three values of the plasmon oscillation permittivity that correspond to 
functions (1271) . are of the form (a = 1, 2, 3) 
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(6+a) 



1 + 



CL\CL2CL'i 



8 





' ' (l-2^)a^(f <i + J T ) 



+4 + s^) £ 

0=1 




-1 



(30) 



Expressions for the electric potential outside a nanoellipsoid for n = 4 can be found from 
(HUD by using Q25H3QD- 

In the case of the higher-order plasmonic modes ( n > 4) the plasmon potentials and the 
resonant permittivity values of a triaxial nanoellipsoid can be found in analogous way. 

To conclude this section, we repeat a general approach to solving a problem of plasmon 
oscillations in the triaxial nanoellipsoid. First, we find polynomial solutions of the Laplace 
equation (the Niven's functions of the respective order). This can always be done. The 
only problem is to find roots of polynomials of several variables. Then using the Dyson's 
theorem we find the integral in Eq.(jHJ) for the eigenvalue. As the integral is known to 
be proportional to the Niven's function it is sufficient to find only the highest terms over 
coordinates, or even one of these terms. The calculation of other terms can be used to 
check calculations only. Then, by finding the ratio of one the highest terms of Niven's 
function to the analogous highest term of the integral we obtain the desired value of the 
resonant permittivity. Using a dispersion law for a particular substance one can find plasmon 
frequencies of the modes under consideration. 

IV. DISCUSSION AND ANALYSIS OF THE RESULTS OBTAINED 

The explicit expressions for plasmon potentials and eigenvalues of the permittivity ob- 
tained in the previous section are of great practical importance. From a general solution for a 
particle in arbitrary external field (j4j) it is seen that the presence of plasmonic nanoparticles 
is manifested in two ways. On the one hand, the spatial field distribution near nanoparticles 
is defined by the eigenfunctions of a corresponding resonant plasmon. So, it is especially im- 
portant to know spatial properties of the plasmon potentials in the region near nanoparticles 



11 



(amplification of local fields for SERS and analogous applications), and far from nanoparti- 
cles (to couple the near- and far- fields, for the nanoantennas and nanolenses). 

The far-field plasmons are characterized by usual multipole moments, so the ALL multi- 
pole moments of ALL plasmonic modes should be known. They may be found from expansion 
of plasmon potentials far from a particle (the mode number is omitted) 



1 3 1 3 1 3 

g (r -> oo) = — ^ d a x a + Q<*P X <* X P + ^7 °c*M x <x x I3 x i + ••• ( 31 ) 

a=l a, (3=1 a,f3,y=l 

where d a , Q a p and O a p~, are the dipole, quadrupole, and octupole moments, respectively 
Multipole moments of different plasmon modes may be found easily by using a operator 



relation between the interior and exterior Niven's functions 



57 



59] 



/ d d d \ J D 2 D 4 11 

9n ~ n {dx~ 1 , dx~ 2 , dx~ 3 ) I 1 + 2 (2ra + 3) + 8 (2n + 3) (2n + 5) + "J r ' (32) 

where H n (x\, x%) is the polynomial derived from the interior Niven's functions using the 

3 

highest coordinate powers only, and where D 2 = ^ a 2 ''" 



a=l 



At large distance from the origin, every differentiation increases the singularity of 1/r. 
The 2 n -th multipole moment of the n-th mode will, therefore, be determined by the principal 
term of Eq.(|32|) 

( 9 9 d \ 1 _ {-l) n {2n)\ 1 
Hn \dx~l' dx~ 2 'dx~Jr- 2 «n! r 2 ^ Hn ^ Xs} ' [dd) 

By comparing that expression with the expansion ( |3TI) one may find the desired 2 n -th mul- 
tipole moments of n-th mode. The multipole moments of higher powers (2 n+1 , 2 n+2 ...) can 
be found from ( 1321) by keeping in series the next terms D 2 , D 4 , . . . The multipole moments 
of lower powers (2 n_1 , 2 n ~ 2 ...) for n-th mode will be equal to zero. 

On the other hand, the amplitude of plasmon resonance is defined by efficiency of a 
plasmon excitation. In its turn this efficiency is product of 2 factors. First is the integral 

of spatial overlapping of external and plasmon field (factor J dr (e n ■ E°) in Eq.(j4])). Second 

v 

factor is determined by a coincidence of external field frequencies with the resonant plasmon 
frequency (factor (e n — e (uj))" 1 in (j4j)). The last one depends considerably on properties of 
a nanoparticle material. The described features of plasmon resonances in nanoparticles will 
be analyzed below. 
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The case n = 1 "dipole" modes 

The electric field of "dipole" modes inside an ellipsoid nanoparticle is homogeneous (fill) . 

and these modes will be most effectively excited by the homogeneous external fields which 

maximaze overlapping integral J dr (e n • E°). In the case of a highly nonhomogeneous ex- 

v 

ternal field, the spacial variations of the field will be averaged and equivalent to a weaker 
homogeneous field. 

A spatial distribution of the squared electric field for modes with n = 1 is shown in 
Fig. 2 where the ellipsoid axes are 02/01 = 0.6 and a 3 /a\ = 0.4. It is seen that for a given 
geometry there are always two maxima localized at a nanoparticle surface. The mode whose 
field is along the longest axis has the highest nonhomogeneity (Fig.2, m — 1). The maximal 
amplification of the external field by that mode is determined by factor 



E 



E° 



1 )e(u) 



(34) 



Figure 3 illustrates the surface charge distribution for modes with n — 1. One can see 
a quasidipole character of those modes. But unlike the spherical plasmons, the plasmons 
of ellipsoid nanoparticles have the high-order multipole moments. For the case n — 1, the 
quadrupole moment of the mode proves to be zero. For dipole and octupole moments we 
obtain (see Eq.( f3TI) ). 

A) For the mode (for more information see Appendix B) 



d a 



B) For the mode g\ 



(2) 



O af3l = <W<V^ - l) I H - 5«a + J] a 2 \ V. 



(35) 



do 



9 



O aPl = SapS^— [ef ] -1)1 2a 2 2 -5c£ + $>2j V. 



(36) 



C) For the mode g{ 



(3) 
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d a = -<5 a3 ^ (ef - l) V, 
Oap, = S aP S 73 ^- (4 3) - l) |2a 2 3 - 5a 2 + £ a\ J V. (37) 

In the case of a sphere (ai = a-i = 03), the dipole mode has a zero octupole moment, as 
follows from (|35l)-(l37|). 

The resonant values of dielectric constant are connected with a corresponding plasmon 
frequency. Above we found the values of resonance permittivity that are independent of a 
nanoparticle material. To calculate the plasmon frequencies it is necessary to use the corre- 
sponding dispersion laws for real substances. Figure 4 demonstrates the plasmon frequency 
dependences of the modes with n = 1 for the Drude's dispersion law u n = Upi/y/1 — e n . It 
is seen that the plasmon frequencies of all modes have a simple monotonic character. As a 
whole, the properties of modes with n = 1 are quite transparent and analogous to dipole 
modes of a sphere. 

The case n = 2 "qudrupole" modes 

Unlike the case with n — 1, the modes with n = 2 may be of two types. They are even 
modes (relative to change of sign of one of the coordinates, two modes) (|T5|) , and odd modes 
( TIB]) (three modes). Such modes cannot be excited by a homogeneous field of any frequency, 
because of a overlap integral is equal to zero for such modes. On the other hand, such modes 
may be excited by fields that are highly nonhomogeneous in the vicinity of a nanoparticle. 

A spatial distribution of a squared electric field of the modes with n = 2 is shown 
in Fig. 5 for an ellipsoid with axes 02/01 = 0.6 and 03/01 = 0.4. From this Figure it is 
seen that for a given geometry there are four maxima and one local minimum of the field. 
Figure 6 illustrates the surface charge distribution for the same modes. It shows vividly a 
quasiquadrupole character of these modes. 

Below we draw down explicit equations for the dipole, quadrupole, and octopole moments 
of the modes with n = 2 (see Eq.( l3Tl) ). The dipole and octopole moments are zero in this 
case. For a quadrupole moment we have the following expressions: 

A) For the mode g% (for more information see Appendix B) 
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,7=1 



B) For the mode 



(2) 



C) For the mode 



v7=l 



g Q/3 = (4 3) - 1) (^2^3^ + 5 a3 <Ws) v. (40) 

D) For the mode 

Q a p = -y^ (4 4) - l) {SaiS&al + 5 a3 ^i«s) V. (41) 

E) For the mode g^ 

Q a p = (4 5) - l) {SaiSpal + <f« 2 fyia*) V. (42) 

Figure 7 demonstrates the dependences of plasmon frequencies of the modes with n = 2 
for the Drude's dispersion law on semi-axis ratio. The plasmon frequencies of the modes 
with n = 2 have a nontrivial character with local maxima and minima. In the case of na- 
noellipsoids of arbitrary shape, the surface ( a 3/ a i ; 02/01) has the form of valleys and 
ridges and has no extremums for nontrivial ellipsoids. From quantum condition that the 
plasmon frequency is proportional to the plasmon energy it becomes clear, from thermody- 
namic viewpoint, that nanoparticles with an excited mode are under strain that is minimal 
at certain geometries. This circumstance could be used in a synthesis of nanoparticles which 
are optimal for the given plasmons. 

From the analysis of Fig. 7 it also follows that under certain ratios of 03/01 (02/01 is fixed), 
the values of plasmon frequencies corresponding to different indices m = 1,2,. ..,5 of the 
multipolarity n = 2, may coincide. So, by using a nonhomogeneous source of electromagnetic 
field one may simultaneously excite several different types (modes) of the plasmon oscillations 
at one frequency. This fact can be used for development of efficient souces of light and related 



devices 



501 ] . When ellipsoid is degenerated into a spheroid (ai = a 2 7^ o 3 ), some of the curves 
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represented in Fig. 7 should coincide. In case of multipolarity n = 2, there should remain 
only n + 1 = 3 independents dispersion curves (and not 2n + 1 = 5, as for an ellipsoid of 
general shape). In a particular case of a sphere ( ) there will remain only one 

dispersion point to be determined by the well known expression w = lu p i a/2/5. 

As a whole, the modes with n — 2, though having much in common with quadrupole 
modes in a sphere, have a variety of nontrivial properties and may be successfully applied 
in different nanooptical devices. 

The case n = 3 "octupole" modes 

The modes with n = 3 are described by Eqs. (l2~2"j) (six modes) and ( |23l) (one mode). They 
can not be excited by a homogeneous field, because one can show that their dipole moments 
are equal to zero. However, such modes may be excited by fields which are nonhomogeneous 
in the vicinity of a particle. 

The squared electric field distribution of these modes is shown in Fig. 8. It is seen that 
the field structure is very complex and has several maxima and minima. Note also that the 
squared field maxima are always on a surface of a nanoparticle. 

A surface charge distribution of modes with n = 3, as shown in Fig. 9, is rather complex 
too. Nevertheless a spatial field structure of those modes can be clearly understood from a 
simultaneous analysis of Figs. 8 and 9. 

Figure 10 illustrates dependences of the plasmon frequencies on one of the ellipsoid semi- 
axes. These dependences also have the maxima and the minima that can be used in a 
number of spectroscopic applications. 

The case n = 4 "hexadecapole" modes 

Spatial structure of the modes with n = 4 is more complex, and we do not consider it in 
this work. We restrict ourselves to dependences of the relative plasmon frequencies / u p i 
(m = 1, 2, 9) on one of the semi-axes 03/01 ration at the given ratio 02/01 = 0.6. (Fig. 11). 
These dependences are similar, as a whole, to the lower-order mode dependences, but now 
one dipsersion curve may have several maxima and minima. 

The case n > 4 
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As it was mentioned already higher orders plasmonic modes in 3-axial nanoellipsoid can be 
found analytically within our general approach. Such analytical solution is very complicated 
and we have restricted our analytical calculations to modes with n < 4 here. 

Nevertheless for the sake of completeness it is important to know (at least qualitatively) 
properties of all higher modes. To calculate higher orders modes we made use of boundary 



elements method (BEM, see e.g. 64J). This method is based on surface integral equation 



2na (Q) = ^ / dS M a (M) ' ^ , (43) 

£ n + 1 Js r MQ 

which can be derived from volume integral equation ([5]). Here a (Q) is surface charge density 
of plasmon mode, rig is normal to surface of nanoparticles at point Q , tmq is radius- 
vector from point M to point Q and integration is over surface of nanoparticle.Then one 
can partition surface 5* of nanoparticle into iV small triangle pieces AS* and reduce surface 
integral equation to system of N linear equations. 

The results of numerical calculations with 3000 triangles are shown in Fig. 12, where 
dispersion curves of 300 modes are shown by circles while solid lines are analytical dispersion 
curves from Figures 4, 7 and 10. 

From this Figure the general structure of dispersion curves of plasmonic modes in triaxial 
nanoellipsoid is clearly seen. In the region 03/01 > 0.1 there is a good agreement between 
analytical and numerical calculations. On the contrary, in the region 0.3/0,1 < 0.1 numerical 
solution is inaccurate and differs substantially from analytical one because of very sharp 
edges of ellipsoid in this region. The analytical solution show that resonat permittivity goes 
to infinity or to zero in the limit of very thin ellipsoid (0.3/01 — > 0). 



V. CONCLUSION 



In this work, the plasmon oscillations of a triaxial nanoellipsoid have been studied in 
a quasistatic approximation. A general approach has been proposed to finding analytical 
expressions for the plasmon electric potentials and the corresponding values of permittivity 
(or plasmon frequencies) of an arbitrary multipolarity order. By making use of that approach 
one could derive the explicit analytical expressions for potentials of several plasmon modes 
and the dielectric permittivity. 
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The analytical expressions obtained in this work are extremely important for nanoplas- 
monics and nanooptics. They can be used to interpret the experimental results on optical 
properties of nanoparticles and to assess reliability of numerical calculations on the interac- 
tion of electromagnetic radiation with nanoparticles that may be approximated by triaxial 
nanoellipsoids. 

The variety and complexity of optical properties of the triaxial nanoellipsoids can also be 
used in many new nanooptical applications. In particular, by choosing a suitable geometry of 
an ellipsoid one can ensure its effective interaction with elementary quantum emitters both 
in the radiation and absorption regimes. As a result, the efficiency of natural and artificial 



fluorophores and other nano dimensional light sources might be substantially enhanced 



50, 



62|. 

In addition, due to the complexity of geometry of an ellipsoid and its plasmon modes, 
one might build, on its basis, the effective optical nanoantennas [lSl^lsj]- They possess all 
properties of more complex in production Yagi-Uda nanoantennas [631 ] . and may be more 
easily synthesized. Of course, a detailed description of a non-dipole radiation pattern of an 
ellipsoid nanoantenna is possible only beyond the framework of quasistatic approximation 
only. However, the results obtained in this work make it possible to perform fast estimations. 
Moreover, within the framework of the longwavelength perturbation theory one can find 
next order corrections (~ 1/A 2 and ~ 1/A 3 ) to resonant values of permittivity and the 
respective eigenfunctions [62]. 
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APPENDIX A: THE EXPLICIT EXPRESSIONS FOR THE INTEGRALS USED 
IN CALCULATIONS OF THE PLASMON POTENTIALS OF A TRIAXIAL NA- 
NOELLIPSOID 



By F.D. Dyson the following theorem has been proved |60j, |6l| (p, q, s — 0, 1, 2, . . .) 



IV I a Is W 2 ! 
,x\ p x 2 q x 3 Tiaras 

\r' — r| 



v 



E 

fc=0 



2 2fc A;! + 1)! 



/ ^— ii- J2 

a=l 



3 2 



fc+1 



i2(u) 



u + ai 



xD 



9 \ P / 1 \ Q / 1 

k i / ajxi \ / 03X2 \ / 03X3 



u + a\) \u + a 2 ./ \' u + a 3, 



(Al) 



where iV=max(p, g, s); -D = ^ " + 2 a " J^- , /)° = 1. Basing on this theorem one can obtain 



a=l 



analytical expressions for the following integrals used in this work (a, (3, 7, a — 1, 2, 3): 



dr' 



r — r 



7raia 2 a 3 / (« ) - a£j a (w ) 



«=i 



(A2) 



r — r 



v 



Tia^a^Xa ( J Q (w ) - ^a^/ a/3 (w ) ) , 

13=1 



(A3) 



|r' — r| 



v 



Tra^a^a^XaXp I I af3 (u ) - ^x^J a/37 (u ) 
V 7=1 

1 / 3 

+-7ia 1 a 2 a 3 a 2 a 5 a(3 I (u ) - 2 a*JJ 7 (u ) 



4 

3 3 

+ x 7 x ^ v Oo) 

7=1 CT = 1 



7=1 



1 / 3 

-naia 2 a 3 a 4 a 5 a(3 I J Q (« Q ) - 2 ^x 7 / Q7 (« Q ) 

V 7=1 



3 3 

+ Yl x j x ° Ia ~t° ( u °"> 

7=1 a=l 



(A4) 
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j ^ \r' — t\ = naia2a 3 a lt a l3 a2 ~l X a X P X l {^aPl ( U o) ~ ^ ^glg^a (^o) 
1 f 3 

+-ira 1 a 2 a 3 5 a/ 3a 2 a a 2 x~ / ( I a (u ) - 2 x 2 a I aa (u ) 

3 3 

£7=1 p=l , 

-^axa 2 a z b a( ia 2 a at f x 1 ( I ay (u ) - 2^2,x 2 I aia (u ) 

3 3 

+ ^2 ^2 x ° x l M 

(7=1 <T=1 y 

+-7raia 2 a 3 5 a/ g5 a7 a^x Q ( J a (u ) - 2^x 2 a l aa (u ) 



(7=1 



(7=1 



2 

3 3 

(7=1 p=l 



(7=1 



1 / 

--7raia 2 a35 a/ 35Q 7 a^x a J aa (u ) - 2^^/ aaCT («o) 



(7=1 



3 3 



+ S S ^P^p H) I , (A5) 

(7=1 p=l / 

in which 5 a p is the Kronecker's delta-symbol equal to unity at a — (3, and equal to zero at 

(3; 



oo 



H u o) = [ ~~ET7~~Xi I a (u )= [ - — TTpTT ' 
J R{u) J (u + a 2 a )R{u) 

UO 

oo 

(«o) = y 

4/37 ( M o) 



MO UO 
OO 
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oo 
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Uo 
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3 

where R (u) = ]~[ y/u + a^; the parameter m = 0, if the point r is inside an ellipsoid, and 

a=l 
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u is the positive root of equation Yl Mo +a 2 = ^ ^ point r is outside the ellipsoid. 



a=l 



APPENDIX B: THE EXPLICIT EXPRESSIONS FOR PLASMON WAVEFUNC- 
TION OUTSIDE ELLIPSOID (SEE EQ.flTUD) 

A) In the case n = 1 we get three functions (a = 1, 2, 3) 



(a) _ I a (Up) 
9l — x a T 



[Bl) 



B) In the case n = 2 we get five functions. Two of them are of the form 



(2) 
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where 
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The rest three have the form 



1/2 



E 
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(B3) 



(3) -^23 (U ) (4) ^31 («o) 
g 2 = X2X3 ~ , g 2 = X3X1 — 



! 23 



'31 



(5) I%2 (U ) 

g\ = x x x 2 — . 

1 12 



C) In the case n = 3we get seven functions. Six of them are (a = 1, 2, 3) 



(B4) 
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and 



93 = x * [ 



(1 + 25 ai ) a* I 
A' - al 
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The rest seventh function is 



(B7) 



(7) 1 123 (Uo) 

g K 3 = x x x 2 x z - 



'123 



(B8) 



In the expressions above we have used notations I a = I a (0), I a p = I a p (0) and J aj a 7 = 7 a/ g 7 (0) 
(see Eq. (1561) ). 

It should be noted, if we set w = (ellipsoid surface), then the potential functions of 
an exterior problem, represented in Appendix B, turn into the respective functions of an 
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interior problem. 
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List of Figure Captions 



Fig.l Geometry of the problem. 

Fig.2 (Color online) Distribution of the squared electric field logarithm (a.u.) near a 
triaxial nanoellipsoid with semi-axes a 2 /a 1 = 0.6 and a 3 /ai = 0.4 in the plane x 3 = for the 
plasmon modes with n — 1. Index m — 1, 2, 3 denotes upper index of the potential function 
(pi . Transition from red color to blue corresponds to transition from the highest squared 
values of electric field to the least squared values. Black line corresponds to the surface of 
the nanoparticle. 

Fig.3 (Color online) Surface charge distribution (a.u.) of plasmon modes of a triaxial 
nanoellipsoid with n — 1 at 02/01 = 0.6 and a 3 /ai = 0.4. Index m — 1,2,3 denotes upper 
index of the potential function (f)^. Warm colors correspond to the positive charge; cold - 
negative. 

Fig. 4 (Color online) Relative frequency of plasmon oscillations of a triaxial nanoellipsoid 
uj^ /u p i as function of semi-axes ratio 03/01 at 02/01 = 0.6. 

Fig. 5 (Color online) Distribution of the squared electric field logarithm (a.u.) near a 
triaxial nanoellipsoid with 02/01 = 0.6 and 03/01 = 0.4 with semi-axes a 2 /oi = 0.6 and 
03/^1 — 0.4 in the plane x 3 = for plasmonic modes with n — 1. Index m = 1, 2, . . . , 5 
denotes upper index of the potential function 4^ ■ Transition from red color to blue cor- 
responds to transition from the highest squared values of electric field to the least squared 
values. Black line corresponds to the surface of the nanoparticle. 

Fig. 6 (Color online) Surface charge distribution (a.u.) of plasmon modes of a triaxial 
nanoellipsoid with n = 2 at 02/01 = 0.6 and 03/01 = 0.4. Index m = 1,2, ... ,5 denotes 
upper index of the potential function (f)^. Warm colors correspond to the positive charge; 
cold - negative. 

Fig. 7 (Color online) Relative frequency of plasmonic oscillations of the triaxial nanoellip- 
soid /(dpi as the function of semi-axes ratio 03/01 at a 2 /ai = 0.6. 

Fig. 8 (Color online) Distribution of the squared electric field logarithm (a.u.) near a 
triaxial nanoellipsoid with semi-axes 02/01 = 0.6 and 03/01 = 0.4 in the plane x 3 = for 
plasmon modes with n — 3. Index m = 1,2, ... ,7 denotes upper index of the potential 
function 3 " 1 ' ) . Transition from red color to blue corresponds to transition from the highest 
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squared values of electric field to the least squared values. Black line corresponds to the 
surface of the nanoparticle. 

Fig.9 (Color online) Surface charge distribution (a.u.) of plasmon modes of a triaxial 
nanoellipsoid with n = 3 at a 2 /oi — 0.6 and 03/01 = 0.4. Index m = 1,2, ... ,7 denotes 
upper index of the potential function (f)^ . Warm colors correspond to the positive charge; 
cold - negative. 

Fig. 10 (Color online) Relative frequency of the triaxial nanoellipsoid plasmon oscillations 
uj^ /uj p i as function of the semi-axes ratio a 3 /oi at a 2 /ai = 0.6. 

Fig. 11 (Color online) Relative frequency of the triaxial nanoellipsoid plasmon oscillations 
ijf 1 ^ juj v i as function of the semi-axes ratio 03/01 at 02/01 = 0.6. 

Fig. 12 (Color online) Permittivity of the triaxial nanoellipsoid plasmon oscillations as 
function of the semi-axes ratio 03/01 at a 2 /ai = 0.6. Solid lines represent analytical solution 
for n — 1 (blue lines), n = 2 (green lines) and n = 3 (red lines). Circles show results of 
numerical calculations within boundary element method. 
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